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Abstract. - We study the effect of electron-phonon interaction on current and zero- frequency 
shot noise in resonant tunneling through a series triple-quantum-dot coupling to a local phonon 
mode by means of a nonperturbative mapping technique along with the Green function formula- 
tion. By fixing the energy difference between the first two quantum dots to be equal to phonon 
frequency and sweeping the level of the third quantum dot, we find a largely enhanced current 
spectrum due to phonon effect, and in particular we predict current peaks corresponding to 
phonon-absorption and -emission assisted resonant tunneling processes, which shows that this 
system can be acted as a sensitive phonon-signal detector or as a cascade phonon generator. 



Phonon-assisted inelastic tunneling in semiconductor quantum dot (QD) system at low 
temperature has become a focus issue in recent years. [1-10] In particular, a recent experiment 
has measured the nonlinear tunneling through a double-QD (DQD) with the observation of 
spontaneous phonon emission leading to an additional satellite peak in the current spectrum, 
[1] which can be ascribed to an interference effect of the electron-phonon interaction (EPI) 
in a DQD via nonperturbative theoretical analyses. [2-4] This experiment opens a possibility 
of designing DQD as a coherent phonon generator. However, the phonon-assisted peak in 
current spectrum is quite fragile and thus detection of phonon-signal is a difficult task in a 
DQD. [1] 

In this letter, we propose a setup containing a triple QD in series coupled to a common 
local phonon bath and two normal leads, in which energy difference between the first two QDs 
is fixed to be equal to the phonon frequency, i.e., head of the device acts as a phonon emitter 
when there is a nonequilibrium current flowing through as suggested by Fujisawa et al. [1] 
Intuitively, it is imaginable that if the energy of the third QD is tuned, via applying gate 
voltage, to be higher than the second QD with one-phonon-energy (case a in Fig. 1 below), 
the emitted phonon could be re-absorbed by electron to help electron tunneling through QD 
3 resonantly, resulting in a phonon-absorption-assisted enhanced peak in current spectrum. 
That is to say that QD 3 detects the generated phonon. On the other hand, we predict a 
significant enhancement of current provided that the energy of QD 3 is further lower than QD 
2 by one-phonon-energy (case c in Fig. 1), showing that more phonon quanta are generated 
in tunneling process. 
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Fig. 1 - Schematic description of tunneling of a triple-QD in the presence of electron-phonon interac- 
tion. e\ — €2 — huio', a, b, and c denotes 63 — £2 = Swo, 0, and —hujo, respectively. Hl([ir) responds to 
chemical potentials of the lead L(lead R). Fl and Tr stands for the coupling between the lead with 
QD. Ti and T2 is the coupling between QDs. 

Our setup is schematically shown in Fig.[TJ We consider a triple QD (consisting of QD 1, 
QD 2, and QD 3) which are connected via a tunnel barrier. QD 1 and QD 3 are connected to an 
electron reservoir in thermal equilibrium with chemical potentials (1l (source) and /ir (drain), 
respectively, with fiL > (1r. Its Hamiltonian can be written as 

H = Hl + Hr + H cen + Ht : (1) 

where Hl + Hr describe the two leads, H cen represents the central region consisting of a triple 
QD in series, and Ht is the coupling of the dots to the lead, respectively: 

H L +Hr= ^2 e r,kclk C r,k, (2) 
rjeL,R;k 

H cen = Y, eotdUa + t^o^b - Y Aa4<U& f + b) - Tx{d\d 2 + 4^) - T 2 {d\d z + d\d 2 ), (3) 

a a 

H T = Y(VLcl k d 1 +V R dic Rk +K.c), (4) 

k 

where c^ fc (c, ; fc) (r/ = {L,R}) and d^da) denote creation (annihilation) operators for spinless 
electrons with momentum k and energy e^k in the left and right leads, and for spinless electrons 
on the ath (a = {1, 2, 3}) QD, respectively. e a is the energy level of the ath QD, and Ti(T 2 ) 
stands for the interdot hopping between the QD 1 and QD 2 (QD 2 and QD 3). Vl<r) is 
the coupling constant of the central regime with lead L(R). This sort of triple-QD system 
was already realized in experiments about ten years ago, in which the energy levels of every 
QD and dot-dot hoppings can all be adjusted by applying gate voltages. [11] The operator 
b^(b) in Eq. ([3]) creates (destroys) a phonon, and A Q is the interaction constant of the electron 
on the ath dot with phonon. Here, we consider the phonon bath as a single-mode phonon 
with dispersionless energy hujQ, and assume that phonon remains coherence during electron 
tunneling processes, i.e., we take no account of phonon dissipation in the present investigation. 
For simplicity, we also ignore electron-electron interaction in this letter. 

As the system under study contains a many-body problem of electron-phonon scattering 
during the single-electron tunneling process, it can not be solved analytically. Applying the 
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nonperturbative mapping technique suggested by Bonca and Trugman [5], we can transform 
the many-body problem into a multichannel one-body problem by transforming the Hamil- 
tonian Eq. ([1]) in terms of a new set of electron states combined of single-electron states and 
n-phonon states. We define a direct product states with single-electron states and n-phonon 
Fock states as: 

\a,n) = dl&-\Q}, (5) 



and 

| 7? fc,n) = ct fe &|0). (<i) 



such that the electron state \a) in QD and the state of electron in leads \rjk) are accompanied by 
n phonons (|0) is the vacuum state). We perform this transformation on the electron-phonon 
interaction part of Eq. ([3]) , so the many body problem can be mapped onto a one-body model: 

-\ a Vn + n + l)(a,n\ + \a,n){a,n + 1|). (7) 

an 

With defining these Dirac bracket as new operators: 

d ln = \ a > n ), c lkn = \vk,Tl), (8) 

we obtain: 

)• (9) 

ctn 

We can do the same transformation on the Hamiltonian Eq.[T]as H = Hi+Hjf+H cen + H-T- [4] 

H L + H R = er,kn4i kn Crikn, (10) 

n£L,R,kn 



an an 



H cen = ^e an dl n d an -^X a Vn + l(d t an+1 d an +}I.c.)-T 1 (d\ rl d 2n +R.c.)-T 2 (dl n d 3n +}I.c.), 

(11) 

H T = ^2(V Ln c[ kn d ln + V Rn d\ n c Rkn + H.c.), (12) 



kn 



where d^ an and c' kn are new operators with phonon quanta n, e an — e a + nfiWQ, e v k n = 
e, q k + nhiOQ. It is important to note that the channel indices stand for the phonon quanta, so 
we have to add a weight factor P n = (1 - e -^o/k B T^ e -nt^ /k B T to tne ntll channel. 

We define the retarded GFs of the triple QDs, G r af}mn {t,t') = -i6(t-t')({d am (t), d\ jn (t')}) 
(a, (3 = 1,2,3), G^ ka mn (t,t') = ~i6(t - t')({c, lkm (t), dl n (t')}) where m,n represents the 
phonon quanta. Employing the equation-of-motion technique, we get the retarded GFs in the 
matrix form: 

<G\ X G\ 2 G[ 3 \ f(u, + iT L )I-A 1 T X I X" 1 

Gil G r 22 G r 23 \ = \ T t I uI-A 2 T 2 I , (13) 

v G 3i G r 32 GlJ \ T 2 I {u: + ^ R )I-A 3 J 

in which to responds to energy, and I is a N x N unit matrix and A a is a N x N symmetrical tri- 
diagonal matrix with A a . nn — e a + nfrui , A a . r ^ n _ l ) = —X a y/n : and A a . n r n+ x) = —X a y/rT+T, 
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respectively. T n — 2nJ2k l^nl ~~ e »?fc — m fi^o) with the wide band limit represents the 
coupling strength of the center region with lead rj. 

The current of the nth channel in the left lead through the center region can be ob- 
tained from the time evolution of the occupation number operator of the left lead: I = 
-e(N L ) = -^([H,N L ]) with N L = Y^k,n c Lkn c Lkn- After some algebra, we find I = 

f Efc,n[^fi{4bi^in) - V Ln( d in c Lkn)}- Using the Keldysh nonequilibrium Green's function 
(GF) technique and considering that the total current is a sum over all pseudo-channels ac- 
companied with the weight factor P n , we get the current as: [4-6] 

1 = { I ^£(4 m nW){P«/2H[l - " Pm/?M[1 - /!»]}, (W) 

ran 

where /l(A)( w ) = (l J t-eS u ~^ Ll - R '>^ kBT )~ 1 is the Fermi distribution of the leads at local thermal 
equilibrium, inwhich T is the temperature. i_L nm represents the transmissionprobability of a 
electron through the center region from the nth channel of the left lead to the mth channel of 
the right lead. According to the Fish-Lee relation relating the scatting matrix elements with 
the retarded GFs, [12] we describe the transmission and reflection probabilities in terms of 
retarded GFs Eq. Ip]): 

rL(R),mn = -S mn + irG[ 1(33) mn (u) , (15) 

tL(R),mn = * rG 31(13),mnM- (16) 

Furthermore, we use the Buttiker scatting method to calculate the shot noise [13] and follow 
our previous work, [4] we obtain the zero-frequency shot noise of the system: 

Sll(0) = ^-[ dcuJ2{\(tW) mn \ 2 P n f 2(^(1 - + \(Ar L ) mn \ 2 [Pnf2(u) 

mn 

x(l - /£(«)) +P m fS(u)(l - fUw))} + \(t R t R ) mn \ 2 P n f^)(l - /j?M)}.(17) 

In order to obtain a more physical view of the results, we now make a numerical simulation 
of the current and the shot noise based on Eqs. ([14"]) and (fl7|) as functions of the energy 
difference between QD 2 and QD 3: e = €2 — £3 (which can be tuned by applying gate voltage). 
With considering the experiment by Fujisawa, [1] we choose the temperature T = 23mk 
{ksT — 2ueV), and dispersionless phonon fkuo = 20ueV. The scale of the dot is d — c s /uo = 
Ylhnm where c s = 5300m/s is the longitudinal sound velocity. At low temperature the 
phonons are assumed to be piezoelectric acoustical mode with interaction constant |A| 2 = 
5tF 2pVu a ' wnere P is the piezoelectric coupling, p is the ion mass density, V is the volume. [14] 
Employing the typical GaAs parameters, [14] we get A ~ 1.6/ieV. We find the electron-phonon 
interaction constant Ai,A2,A3 on the three quantum dots coincide up to a phase factor. [2] 
For simplicity, we choose phase factor is —1 between the QD 1 and QD 2 (QD 2 and QD 
3), so we get Ai = — A 2 = A 3 = 1.6ucV (This can be done by tuning size of each dot and 
distance between them in experiments). In this cases, we set the couplings of the dots with 
the two leads to be symmetric, Tl=T r = T~ AneV, and the interdot hoppings to be equal, 
Ti = T2 = T c = AfieV. We assume e 2 = 0/ieV, which is equal to the Fermi energies of the 
left lead and right lead at equilibrium. We also set the energy gap of the QD 1 and QD 2 as 
ei — e 2 = hu>o = 20[ieV and the symmetrically applied bias voltage (J.l — n R — eVsd = 200/ieV 
so that /i£ 3> ei, £2, £3 3* ^r, in which condition sufficient more pseudo-channels are involved 
in transport and thus our numerical results are nearly independent of the bias-voltage. 
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Fig. 2 - Calculated total tunneling current versus energy difference e with parameter: F = 4/ieV and 
Ai = — A2 = A3 = 1.6/ieV, the temperature of the system T = 23mK. 



In Fig[2]we plot the total current of the triple-QD system, for comparison, in the absence of 
EPI (dash line) vs the energy difference e. It is interesting to find a peak at e = e 2 — £3 = — fkvo 
(where ei = £3) in the system without considering electron-phonon interaction. This peak is 
due to a resonant tunneling between QD1 and QD3 via QD 2. We can describe the dynamics of 
the triple QD system with quantum rate equations for the density matrix of the dots. [15] With 
solving the master equation algebraically, we get the current at zero temperature analytically 
as: 

e 4TT 4 
I= h 16T4 + (2r 2 + 8£ 2 + 16££ 12 + 8e 2 12 )T? + 4ef 2 + e\ 2 T 2 + 4e 2 2 e 2 + 8e? 2 e' (18) 

with e 12 = e% — e 2 . It is obvious to obtain the peak of current at e = — e 12 = Hloq. 

In Figj2] we also plot the total current of the triple-QD system in the presence of EPI vs 
the energy difference e. Similar with the current spectrum of a double-QD system, [1-4] an 
obvious overall enhancement of current is found for the triple-QD system with phonon bath 
in comparison with that without EPI. Moreover, as above expected, our results predict three 
current peaks labeled by a, b, and c, which are corresponding to the three configurations 
described in Fig. 1. When an extremely large bias- voltage is applied to the whole system, 
an electron can resonantly jump from QD 1 to QD 2 accompanied by one-phonon generated, 
then if the system is at configuration b, e 2 = 63, the electron can directly resonantly tunnel 
through the device and the current peak b is caused by phonon emission between QD1 and 
QD2; however, if the system is at configuration c, e 2 — €3 = fvjjQ, the electron will emit 
another phonon to reach QD 3, leading to two-phonon-emission peak c; more interestingly, 
if the system is at configuration a, e 2 — £3 = —fkoo, an enhanced current peak a also occurs 
stemming from re-absorption of the generated phonon in head of the device by electron to 
overcome the barrier between QD 2 and QD 3, which is enhanced by 30% in comparison to 
the resonant peak without EPI. Of course, phonon dissipation to environment could lessen 
these phonon-assisted resonant peaks. Albeit we do not consider nonequilibrium phonon 
effect [7,8] and dissipation due to environment, [9,10] however, it can be expected from our 
above calculations that the multi-QD system could function as a cascade phonon generator 
with higher emission efficiency than a double-QD or as a sensitive and "good" phonon-signal 
detector. Actually, we can put this system into a phonon resonant cavity in experiments to 
lessen the leakage of phonon and thus to enhance the probability of re-absorption of phonon. 
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Fig. 3 - (a) Zero-frequency shot noise versus energy difference e with the same parameter of Fig[2] 
(b) Fano factor versus energy difference e. 



These studies will be the subject in our future work. 

We also examine the zero-frequency shot-noise (a) and Fano factor F — Sll/^I (b) with 
the energy difference e, as plotted in Fig|3j Corresponding to the three current peaks due 
to phonon-absorption or -emission, the shot noise spectrum shows also single peak in the 
phonon- assisted resonant tunneling region. As the shot noise S ~ Tt r (l — Tt r ), which T± r is 
transmission probability, there is no shot noise generation when the T tr = 1 or 0. The maximal 
generation of the shot noise occurs while the transmission probability is between and 1. We 
also find a small peak at e = in the absence of phonon. However, the shoulder at e = is 
smeared smoothly by the interaction of phonon. The Fano factor F displays three dips at the 
three points a, b, and c, implying that the inelastic resonance suppresses the shot noise. It 
is worth noticing that the phonon-absorption-assisted tunneling induces a more pronounced 
dip. 

In summary, we investigate the resonant tunneling through a triple QD in the presence of 
EPI at low temperature by means of the nonperturbative mapping technique in combination 
with the nonequilibrium GF method. By making the first two dots act as a phonon emitter 
and sweeping the energy level of the third QD, we evaluate the tunneling current and predict 
that resonant peaks in the current spectrum are not only due to spontaneous phonon-emission 
assisted tunneling, but also due to phonon-absorption process. We also study zero-frequency 
shot noise of this sytem. 

This work was supported by Projects of the National Science Foundation of China, the 
Shanghai Municipal Commission of Science and Technology, the Shanghai Pujiang Program, 
and NCET. 
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